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ABSTRACT 

", The emission of 21-cm radiation from a neutral intergalactic medium (IGM) at high 

ON I redshift is discussed in connection with the thermal and ionization history of the uni- 

verse. The physical mechanisms that make such radiation detectable against the cosmic 
microwave background include Lya coupling of the hydrogen spin temperature to the 
kinetic temperature of the gas and preheating of the IGM by the first generation of stars 
\q . and quasars. Three different signatures are investigated in detail: (a) the fluctuations in 

the redshifted 21-cm emission induced by the gas density inhomogeneities that develop 
| at early times in cold dark matter (CDM) dominated cosmologies; (b) the sharp absorp- 

tion feature in the radio sky due to the rapid rise of the Lya continuum background 
that marks the birth of the first UV sources in the universe; and (c) the 21-cm emission 
and absorption shells that are generated on several Mpc scales around the first bright 
quasars. Future radio observations with projected facilities like the Giant Metrewave 
Radio Telescope and the Square Kilometer Array may shed light on the power spectrum 
of density fluctuations at z > 5, and map the end of the "dark ages", i.e. the transition 
from the post-recombination universe to one populated with radiation sources. 
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1. Introduction 

The development of structure in the universe was well advanced at high redshifts. Quasars have 
been detected nearly to z = 5, and the most distant galaxies to even greater distances (Weymann et 
al. 1998). The Lya forest in the spectra of high redshift QSOs shows that the intergalactic medium 
(IGM) itself developed extensive nonlinear structures by early times. Still unknown, however, is 
the epoch during which the first generation of stars was formed, and the time scale over which the 
transition from a neutral universe to one that is almost completely reionized took place. While 
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numerical simulations have shown that the IGM is expected to fragment into structures at early 
times in cold dark matter (CDM) cosmogonies (Cen et al. 1994; Zhang et al. 1995; Hernquist et 
al. 1996; Zhang et al. 1998), the same simulations are much less able to predict the efficiency with 
which the first gravitationally collapsed objects lit up the universe at the end of the "dark age" 
(Rees 1998). The nature of the first dominant ionizing sources is still unknown, and although IR 
observations will permit even higher redshift galaxies to be observed, detections become increasingly 
difficult because of the surface brightness dimming from cosmological expansion. 

Madau, Meiksin, & Rees (1997, hereafter MMR) demonstrated that radiation sources prior 
to the epoch of reionization may be detected indirectly through their impact on the surrounding 
neutral IGM and the resulting emission or absorption against the cosmic microwave background 
(CMB) at the frequency corresponding to the redshifted 21-cm line. While its energy density is 
estimated to be about two orders of magnitude lower than the CMB, the 21-cm spectral signature 
will display angular structure as well as structure in redshift space due to inhomogeneities in 
the gas density field (Hogan & Rees 1979; Scott & Rees 1990), hydrogen ionized fraction, and 
spin temperature (MMR). Because of the smoothness of the CMB sky, fluctuations in the 21-cm 
radiation will dominate the CMB fluctuations by about 2 orders of magnitude on several arcmin 
scales. Radio observations at meter wavelengths have therefore the potential to tell us how (as well 
as when) the primordial gas was reheated and reionized: if the reheating was due to ultraluminous 
but sparsely distributed sources like QSOs, or by fainter but more uniformly distributed objects 
like an abundant population of pregalactic stars. 

In this paper we discuss possible ways in which the next generation of radio telescopes could 
provide the first observations of a neutral IGM at a level of a few tenths of fiJy arcmin -2 . The 
proposed Square Kilometer Array (SKA) would provide an enhancement of about two orders of 
magnitude in sensitivity relative to present-day radio telescopes (Taylor & Braun 1999). We discuss 
the technical requirements necessary to detect the 21-cm signature from the IGM and describe a 
possible design for a special purpose telescope. We emphasize that within existing technological 
limits, radio astronomy could effectively open up much of the universe to a direct study of the 
reheating and reionization epochs. 

2. Basic theory 

The ground state hyperfine levels of hydrogen will very quickly reach thermal equilibrium with 
the CMB. In order to detect the IGM in either emission or absorption, it is crucial that an effective 
mechanism be present to decouple the hydrogen spin states from the CMB. Following MMR, we 
briefly review below the mechanisms that determine the level populations of the hyperfine state of 
neutral hydrogen in the intergalactic medium (see also Meiksin 1999). 
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2.1. Spin temperature 

The emission or absorption of 21-cm radiation from a neutral IGM is governed by the hydrogen 
spin temperature Ts, defined by 

^ = 3exp(-T,/T 5 ), (1) 
n 

where no and n\ are the singlet and triplet n = 1 hyperfine levels, T* = hrvio/kp = 0.07 K, 
i^io is the frequency of the 21-cm transition, hp is Planck's constant, and kp is Boltzmann's 
constant. In the presence of only the CMB radiation with Xcmb = 2.73 (1 + z) K, the spin states 
will reach thermal equilibrium with the CMB on a timescale of ^/(Tcmb^io) ~ 3 X 10 5 (1 + yr 
(Aio = 2.9 x 10~ 15 s _1 is the spontaneous decay rate of the hyperfine transition of atomic hydrogen), 
and intergalactic H I will produce neither an absorption nor an emission signature. A mechanism is 
required that decouples Ts and TcmBj e.g. by coupling the spin temperature instead to the kinetic 
temperature Tk of the gas itself. Two mechanisms are available, collisions between hydrogen 
atoms (Purcell & Field 1956) and scattering by Lya photons (Wouthuysen 1952; Field 1958). The 
collision-induced coupling between the spin and kinetic temperatures is dominated by the spin- 
exchange process between the colliding hydrogen atoms. The rate, however, is too small for realistic 
IGM densities at the redshifts of interest, although collisions may be important in dense regions, 
5p/p > 30[(1 + z)/10]~ 2 (MMR). 

In the low density IGM, the dominant mechanism is the scattering of continuum UV photons 
redshifted by the Hubble expansion into local Lya photons. The many scatterings mix the hyperfine 
levels of neutral hydrogen in its ground state via intermediate transitions to the 2p state, the 
Wouthuysen-Field process. An atom initially in the n = 1 singlet state absorbs a Lya photon that 
puts it in an n = 2 state, allowing it to return to the triplet n = 1 state by a spontaneous decay. 
As the neutral IGM is highly opaque to resonant scattering, the shape of the continuum radiation 
spectrum near Lya will follow a Boltzmann distribution with a temperature given by the kinetic 
temperature of the IGM (Field 1959). In this case the spin temperature of neutral hydrogen is a 
weighted mean between the matter and CMB temperatures, 

rp Tqmb + y a TR fn , 

Ts = i + y~ a ' (2) 

where 

y a = (3) 

is the Lya pumping efficiency. Here P\q is the indirect de-excitation rate of the triplet n = 1 
state via absorption of a Lya photon to the n = 2 level, and Ts 3> T* is assumed. If y a is large, 



5 In the presence of a radio source, the antenna temperature of the radio emission should be added to Tomb in eq. 
(S). The radio emission may make an important contribution in the vicinity of a radio-loud quasar (Bahcall & Ekers 
1969), and would itself permit the IGM to be detected in 21-cm radiation. This may be an important process if the 
emission of ionizing radiation is more tightly beamed than the radio continuum emission, allowing the gas to remain 
neutral near the quasar after it has switched on (MMR). 
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Ts — ► signifying equilibrium with the matter. A consideration of the net transition rates 
between the various hyperfine n = 1 and n = 2 levels above shows that the 1 — ► transition rate 
via Lya scattering is related to the total rate P a by Pio = 4P a /27 (Field 1958). This relation and 
equation (^) are derived in the Appendix. In the limit Tk 3> Tqmb, the fractional deviation in a 
steady state of the spin temperature from the temperature of the CMB is 

Ts — Tcmb ^ A | 7cmb \ -1 
Ts \ y a T K ) 

There exists then a critical value of P a which, if greatly exceeded, would drive Ts — > Tk- This 
"thermalization" rate is (MMR) 

27Ai Tcmb . . fn n . , 1 „_12 _-u ( 1 + z 



Pth = M p.6 x «ru _ | . (5) 



2.2. Brightness temperature 

Consider a patch of neutral hydrogen with an angular diameter on the sky larger than the 
beamwidth of the telescope, and a radial velocity width broader than the bandwidth. In an 
Einstein-de Sitter universe with Hq = 100 h kms _1 Mpc _1 and baryon density parameter Qf,, the 
optical depth of the patch at 21(1 + z) cm for a spin temperature Ts ^ Tcmb is 

T(l) = J_ A?oAlo Ik!W2) (1 + 2)1 .^ 0.0033/,-' (*Z2) (1±±)". (6) 

v ) 327r io w Ts Hq \ ) \ T s ) \0.02 J V 10 J 

The optical depth will typically be much less than unity. Since the brightness temperature through 
the IGM is T^ = TcMBe -T + Tg(l — e~ T ), the differential antenna temperature observed at Earth 
between this region and the CMB will be 

5T b = (1 + z)-\T s - T CMB )(1 - O « (9.0 mK)h^ (i±£) ^ (^f^) ■ (7) 

If the intergalactic gas has been significantly preheated, Ts will be much larger than Tcmb> and 
the IGM will be observed in emission at a level that is independent of the exact value of Ts- By 
contrast, when Tcmb 3> Ts (negligible preheating), the gas will appear in absorption, and \STf,\ 
will be a factor ~ Tqmb/Ts larger than in emission, so that it becomes relatively easier to detect 
intergalactic H I (Scott & Rees 1990). 

The presence of a sufficient flux of Lya photons will thus render the neutral IGM "visible." 
Without heating sources, the adiabatic expansion of the universe will lower the kinetic temperature 
of the IGM well below that of the CMB, and the IGM will be detectable through its absorption. If 
there are sources of radiation that preheat the IGM, it may instead be possible to detect the IGM 
in emission. 
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2.3. Preheating 

The energetic demand for heating the IGM above the CMB temperature is meager, only 
~ 0.004 eV per particle at z ~ 10. Consequently, even relatively inefficient heating mechanisms 
may be important warming sources well before the IGM has been reionized.[] Possible preheating 
sources are soft X-rays from an early generation of quasars or thermal bremsstrahlung emission 
produced by the ionized gas in the collapsed halos of young galaxies (MMR) . While photons emitted 
by a quasar just shortward of the photoelectric edge are absorbed at the ionization front, photons 
of much shorter wavelength will be able to propagate into the neutral IGM much further. Most 
of the photoelectric heating of the IGM by a QSO is accomplished by soft X-rays. The timescale 
required for the radiation near the light front to heat the intergalactic gas to a temperature above 
that of the CMB is typically 10% of the Hubble time. The H n region produced by a QSO will 
therefore be preceded by a warming front. 

An additional heating source is the Lya photon scattering itself (MMR). The heating rate by 
the recoil of scattered Lya photons is 

K = (^-) h P u a P a , (8) 

where P a is the Lya scattering rate per hydrogen atom, and mn is the mass of a hydrogen atom. 
In the case of excitation at the thermalization rate Pth> equation (||) becomes 



l + z 
10 



(9) 



The characteristic timescale for heating the medium above the CMB temperature via Lya resonant 
scattering at this rate is 

about 20% of the Hubble time at z = 8. This result has an important consequence: if a substantial 
population of sources of UV radiation turn on quite rapidly at (say) z ~ 10, there exists a finite 
interval of time during which Lya photons couple the spin temperature to the kinetic temperature 
of the IGM before heating the gas above the CMB (reionization occurs at even later times). A 
window is then created in redshift space where a large fraction of intergalactic gas may be observed 
at ~ 160 MHz in absorption against the CMB. We will show next how this transient phase may 
imprint a strong signature on the radio sky. 



3. Radio signatures 



In an inhomogeneous universe reheating appears to be a slow process occurring at early epochs, followed by 
sudden reionization at later times (Gnedin & Ostriker 1997). 
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3.1. The "cosmic web" in emission at 21 cm 

Numerical ^-body/hydrodynamics simulations of structure formation within the framework 
of CDM-dominated cosmologies (Cen et al. 1994; Zhang et al. 1995; Hernquist et al. 1996; Zhang 
et al. 1998) have recently provided a definite picture for the topology of the IGM at z < 5, one of 
an interconnected network of sheets and filaments with virialized systems located at their points 
of intersection. In between the filaments are underdense regions. This "cosmic web" appears to be 
a generic feature of CDM (Pogosyan et al. 1998), a pattern imprinted in the matter fluctuations 
at decoupling and sharpened by gravity. In a UV photoionizing background, modestly overdense 
filaments (1 < Pb/pb < 5) will give rise to the Lya forest seen in the spectra of high redshift quasars. 
As discussed below, radio observations at meter wavelengths on several arcmin scales may probe 
the cosmic web at times prior to the epoch of reionization. 

Because of the dumpiness of the IGM, the 21-cm optical depth of equation (||) will vary for 
different lines of sight. To estimate the rms level of brightness temperature fluctuations at 21(l + z) 
cm, we shall consider three different cosmological models with parameters suggested by a variety 
of recent observations: an open model (OCDM, £Im = 0.4, h = 0.65, n = 1), a flat model (ACDM, 
Qm = 0.3, S7a = 0.7, h = 0.7, n = 1), and an Einstein-de Sitter tilted model (tCDM, O^f = 1, 
S!a = 0, h = 0.5, n = 0.8). In all cases the amplitude of the power spectrum or, equivalently, the 
value of the rms mass fluctuation in a 8 h" 1 Mpc sphere, ag, has been fixed in order to reproduce 
the observed abundance of rich galaxy clusters in the local universe (e.g. Eke et al. 1998). The 
tilted model has been designed to also match COBE measurements of the CMB on large scales. 
Our choice of an open CDM cosmology (with a normalization which is 2 rms away from the value 
observed by COBE) is guided by the growing evidence in favour of a low value of Qm from cluster 
studies (e.g. Carlberg et al. 1997). The A-dominated model matches both the cluster constraint on 
small scales and COBE constraint on large. It is also able to account for most known observational 
constraints, ranging from globular clusters ages to CMB anisotropies (Ostriker & Steinhardt 1995) 
to recent analyses of the Type la SNe Hubble diagram (Perlmutter et al. 1998; Riess et al. 1998), 
and just fits within the constraints imposed by gravitational lensing (Kochanek 1996; Cooray et 
al. 1999). In all cosmologies the baryon density is Q^/i 2 = 0.02 (Buries &: Tytler 1998). The model 
parameters are sumarized in Table ffl. 



Table 1. Cosmological parameters. 



Model Qm h n as 



tCDM 1.0 
OCDM 0.4 
ACDM 0.3 



0.0 0.5 0.8 0.55 
0.0 0.65 1.0 0.92 
0.7 0.7 1.0 1.1 
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As an illustrative example, consider a scenario in which sources of Lya photons are in sufficient 
abundance throughout the universe to preheat the IGM to a temperature well above that of the 
CMB, and to couple the spin temperature to the kinetic temperature of the intergalactic gas 
everywhere. While the average 21-cm signal is two orders of magnitude lower than that of the 
CMB (and will be swamped by the much stronger non-thermal backgrounds that dominate the 
radio sky at meter wavelengths and that must be removed), its fluctuations relative to the mean of 
the surveyed area - induced in this example by density inhomogeneities only - will greatly exceed 
those of the CMB. Brightness temperature fluctuations will be present both in frequency and in 
angle across the sky. At high redshifts, and with a spatial resolution of ~ 1 arcmin, a patchwork 
of radio emission will be produced by regions with overdensities still in the linear regime and by 
voids. From equations (||) and (|7|), the rms temperature fluctuation relative to the mean, neglecting 
angular correlations in the projected IGM density, is 

(5T*) 1 / 2 « (2.9mK)r 1 f^!)- J l + z) ^ p N9 — , (11) 

V b/ V ' V 0.02 /[0 M (l + z)3 + Q^(l + z) 2 + A ] 1 /2' v > 

where fix = 1 — &m ~ is the curvature contribution to the present density parameter. This 
corresponds to a rms fluctuation in the observed flux of 

Here a p is the variance of the density field in a volume corresponding to a given bandwidth 
[Au/u = Az/(1 + z)\ and angular size AO. The former corresponds to a comoving length L ~ 
(1 + z)cH(z)^ 1 Au/u Mpc, where H(z) is the Hubble constant and v = 1-4/(1 + z) GHz is the 
observation frequency, the latter to a comoving radius R = O.5A0(1 + z)/d a (z), where d a is the 
angular diameter distance. 

On the scales of interest here, one can use the linearly evolved power spectrum P(k) (Bardeen 
et al. 1986) and a filter function that is a cylinder of radius R and height L, to compute the density 
variance: 



p tt 2 R 2 L 2 J Jo x 2 (l-x 2 ) y J ' k 2 

where J\ is the Bessel function of order one, D(z) is the linear growth factor (e.g. Peebles 1993), and 
/ ~ Q^f takes into account the effect of the peculiar velocities in compressing the same emitting 
volume into a narrower bandwidth (as the Hubble flow is reduced in overdense regions, see Kaiser 
1987). In this regime it is also legitimate to assume that the baryons follow the dark matter, since 
segregation effects are important only on smaller scales. We have verified using iV-body simulations 
the validity of the linear approximation and the neglect of angular correlations on scales of ~ 1 
arcmin. 

Figure 1 shows the fluctuation levels at 150 MHz (z = 8.5) as a function of beam size for 
the models listed in Table 1. At a fixed bandwidth of 1 MHz, (5T 2 ) 1 / 2 increases with decreasing 
angular scale from about 1 to 10 mK as a p increases with decreasing linear scale. Since the growth 
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of density fluctuations ceases early on in an open universe (and the power spectrum is normalized 
to the abundance of clusters today), the signal at a given AO is much larger in OCDM than in 
tCDM at high redshift. In a ACDM universe the level of fluctuations is comparable with OCDM. 
The range of density fluctuations that would be detectable at the 5cr level by a projected facility 
like the SKA is shown in Figure 2, as a function of beam size and frequency bandwidth, and for 
several assumed integration times. The detection thresholds are scaled according to a theoretical 
rms continuum noise in a 80 MHz bandwidth at 150 MHz of 64 nJy over an 8 hour integration, 
allowing for only 21 of the 30 dishes in a compact array configuration (Taylor & Braun 1999)0, 



The central compact array in the 'straw man' SKA design proposal has a diameter of ~ 50 km in 
order to achieve 1-10 arcsec resolution. A more appropriate resolution for measuring fluctuations 
in the IGM intensity is ~ 1 arcmin. The compact array would have then to be scaled down to 
D c ~ 7 km (A#/arcmin) -1 to be optimized for such an experiment. The dashed lines show curves 
of constant rms fluctuation of the signal in the corresponding beam. Each angular resolution 
observation corresponds to a different compact array diameter D c (in other words, the array is 
being optimally resized to match its resolution to the desired angular scale). Because the received 
flux is proportional to the solid angle of the beam, the signal decreases with decreasing angular 
scale until it falls below the detection threshold, indicated by the solid lines. With 100 hours of 
observation in a tCDM model, 2 fiJy fluctuations will be detectable only on scales exceeding 3 
arcmin. In OCDM, and similarly in ACDM, a signal of ~ 2 /iJy may be detected at significantly 
higher resolution with the same integration time. 

Figure 3 show radio maps at z = 8.5 for tCDM and OCDM, assuming a beamwidth of 2 arcmin 
and a frequency window of 1 MHz around 150 MHz. Here a collisionless A r -body simulation with 
64 3 particles has been performed with Hydra (Couchman, Thomas, & Pearce 1995). The simulation 
box size is 20h~ 1 comoving Mpc, corresponding to 17 (11) arcmin in tCDM (OCDM). The baryons 
are assumed to trace the dark matter distribution without any biasing. From a visual inspection 
of the two simulated images, as well as from Figure 2, it seems possible that observations with a 
modified (much less spread out) SKA design may be used to reconstruct the matter density field 
at redshifts between the epoch probed by galaxy surveys and recombination, on scales as small as 
0.5 — 2 h^ 1 comoving Mpc, i.e. masses in the range between 10 12 and 10 13 1i~ 1 Mq. 



In the real universe, the observed radio pattern on the sky depend on a combination of the 
underlying Gaussian density field and the distribution and luminosities of the sources that cause 




(14) 



3.2. Observing the epoch of the first stars 



' Sec also 



http: / /www. nfra.nl/skai/science/ . 
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the preheating and reionization of the IGM. In this and the following sections we shall discuss the 
detectability of the first radiation sources in the universe through their impact on the surrounding 
neutral IGM. 

In hierarchical clustering theories for the origin of cosmic structure, massive objects grow 
'bottom-up' - low mass perturbations collapse early, then merge into progressively larger systems. 
In such a scenario, one naturally expects the high redshift universe to contain small scale substruc- 
tures, systems with a virial temperature well below 10 4 K. The resultant virialized systems will 
remain as neutral gas clouds unless they can cool due to molecular hydrogen, in which case they 
may form Pop III stars. These stars will produce a background of UV continuum photons near the 
Lya frequency which will immediately escape into the IGM, while at the same time ionizing the 
gas in their vicinity. Lya photons will propagate into uncollapsed, largely neutral regions of the 
IGM where the kinetic temperature in the absence of preheating will be 

T K {z) « 26 mK (1 + zf (15) 

(Couchman 1985), well below Tcmb because of adiabatic cooling during cosmic expansion. If Pop 
III sources are in sufficient abundance throughout the universe, the Lya flux will couple the spin 
temperature to the kinetic temperature, and T$ will be pulled below Tcmb everywhere. After 



about 10 yr (see eq. 10), the same Lya photons responsible for the coupling will warm the IGM 
to a temperature well above that of the CMB. Eventually a sufficient number of photons above 1 
Ryd may be able to penetrate all of the surrounding neutral gas in the cloud, reach the external 
intergalactic gas, and create expanding intergalactic H n regions. 

The crucial assumption we make here is that the Lya radiation field switches on at a redshift 
z > zth following the formation of the first stars, and reaches the thermalization value Tth at z th 
on a timescale much shorter than the Hubble time at that epoch. While in hierarchical models of 
structure formation the spatial number density of halos of a given mass grows initially exponentially 
(Press & Schechter 1974), the rise time of the UV radiation background will be modulated by 
gasdynamical processes. In the model of Gnedin & Ostriker (1997), for instance, a Lya metagalactic 
flux is rapidly produced in the range 7 < z < 10. In this case coupling between T$ and Tr <C Tcmb 
is almost instantaneous at z = z t h> and the IGM is suddenly detectable in absorption. For z < z^, 
Lya photons begin to heat the IGM and the kinetic temperature increases at the rate 

dT K _ 2n E th dt ^ 2 T K 
dz 3 fee dz (1 + z) 

until Tk > ?cmb and the IGM is now detectable in emission. Here, fi = 16/13 is the mean 
molecular weight for a neutral gas with a fractional abundance by mass of hydrogen equal to 0.75. 
Such a chain of events will leave a strong imprint on the CMB, and mark the epoch of reheating by 
the first generation of stars. At 150 MHz (z t h = 8.5), the brightness temperature (shown in Figure 
4) shows an absorption feature of 40 mK, with a width of about 10 MHz. Such a depression in 
the CMB is much easier to observe compared with brightness fluctuations, as the isotropic signal 
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may be recorded with a larger beam area. Assuming AT ~ 40 mK as a fiducial estimate, the 
signal-to-noise ratio for the SKA is 

N ~ vlMrfeJ \ 100 hours/ b ' { ' 

where Au is the frequency resolution, At, is the beam area in arcmin 2 , and r is the observation time. 
While the strength of this feature depends mainly on the timescale i t h ° ver which the continuum 
Lya background field reaches the thermalization value, we have checked that the absorption dip 
is well above the detection limit for all timescales tth ^ 0.1i#, where tu is the then Hubble time. 
Note that similar sharp signatures in the radio background at meter wavelengths could also be 
produced at the reionization epoch, as recently discussed by Gnedin & Ostriker (1997) and Shaver 
et al. (1999). 



3.3. Observing the epoch of the first quasars 

Very little is actually known observationally about the nature of the first bound objects and 
the thermal state of the universe at early epochs. If ionizing sources are uniformly distributed, like 
an abundant population of pregalactic stars, the ionization and thermal state of the IGM will be 
the same everywhere at any given epoch, with the neutral fraction decreasing rapidly with cosmic 
time. As discussed above, this may be the case in CDM dominated cosmologies, since bound objects 
sufficiently massive (~ 1O 6 M0) to make stars form at high redshift. On the other hand, reionization 
may also occur in a highly inhomogeneous fashion, as widely separated but very luminous sources 
of photoionizing radiation such as QSOs, present at the time the IGM is largely neutral, generate 
expanding H n regions on Mpc scales: the universe will be divided into an ionized phase whose 
filling factor increases with time, and an ever shrinking neutral phase. If the ionizing sources are 
randomly distributed, the H n regions will be spatially isolated at early epochs. 

In such a scenario, 21-cm emission on Mpc scales will be produced in the quasar neighborhood 
as the medium surrounding it is heated by soft X-rays from the QSO itself (MMR). Figure [| shows 
the (differential) radio map resulting from a QSO 'sphere of influence' soon after it turns on at z = 
8.5 (tCDM). The kinetic temperature profile and radiation field around the quasar were computed 
including finite light travel time effects and heating by secondary electrons collisionally produced 
by primary electrons photo-ejected by the soft X-rays (MMR). The specific luminosity of the QSO 
was taken to be proportional to z^~ L8 , normalized to a production rate of H I ionizing photons of 
10 57 s . The IGM was assumed uniform with = 0.08 and Hq = 50kms _1 Mpc -1 . The resulting 
radial temperature profiles were then superimposed on the surrounding density fluctuations as 
computed using Hydra. The visual effect is due to the convolution of the spin temperature profile 
with the (linearly) perturbed density field around the quasar. Outside the H n bubble, there is 
an inner thin shell of neutral gas where the IGM is heated to T$ = Tk > Tqmb as the hyperfine 
levels are mixed by Lya continuum photons from the QSO, and a much larger external shell 
where T$ = Tjc < Tcmb because of adiabatic cooling. At larger distances from the quasar the 
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Lya coupling strength is weakening as r~ 2 , and T5 — ► Tcmb- As the warming front produced by 
the quasar expands, a growing amount of the surrounding IGM is unveiled both in emission and in 
absorption. The size and intensity of the detectable 21-cm region depend on the QSO luminosity 
and age. In the figure the signal ranges from — 3/iJy to 3[iJy per beam. 

Note that, although the quasar was placed in a high density region in the corner of the simu- 
lation volume, the figure can equally be viewed as the emission due to heating by a beam of soft 
X-rays from the quasar with an opening angle of 90°. Thus imaging the gas surrounding a quasar 
in 21-cm emission could provide a direct means of measuring the opening angle of quasar emission. 

4. A compact, special purpose, low frequency array 

The requirements for observing the epoch of IGM reheating and reionization are a narrow 
bandwidth (Af « 1 MHz), to achieve adequate resolution in velocity space at the relevant redshift, 
and an angular resolution of a few arcminutes, the typical scale of the expected density fluctuations. 
The straw man proposal for the SKA suggests a large diameter for its central compact array in order 
to achieve a resolution of ~ 1 — 10 arcsec. While gas density fluctuations would persist to smaller 
scales in the IGM, Figure 2 shows that the associated signal would drop well below the threshold 
for a practical detection. For a successful experiment it is necessary to significantly decrease the 
diameter of the compact array. The current design for the SKA (and the GMRT) has too widely 
distributed 'sub-apertures' for a program oriented toward detecting arcmin scale structures. 

An alternative possibility is to build a compact, special purpose radio telescope dedicated to 
searching for the expected reheating and reionization signals at low frequencies. A single dish with 
a diameter of 200 m (smaller than Arecibo) would resolve scales of 30 arcminutes, corresponding 
to 35 — 50/i _1 comoving Mpc. This is comparable to the largest scale expected for coherent pri- 
mordial density fluctuations (Pogosyan et al. 1998). Figure 1 shows that the expected brightness 
temperature fluctuations on this scale will be in excess of 1 mK. The rms brightness temperature 
is given by the radiometer equation as 

( ^svs A /lMHzX 1 / 2 / 100 hours \ ^ 2 , , 

Aft - (0.33 m K)(^)(— ) (__) (18) 

(Burke & Graham-Smith 1997) for a bandwidth Av and integration time r. The system tempera- 
ture T sys is sky-limited to 200 K in the coldest directions through the Galaxy (Burke & Graham- 
Smith) . If a density fluctuation filled the beam, it would be detectable at the 5a level in a 300 hour 
integration. 

Since such structures tend to be filamentary, a more typical covering factor may be ~ 0.3, 
increasing the required integration time by a factor of 10. This may be offset by increasing the 
bandwidth. The velocity width of a 50/i _1 comoving Mpc filament at z = 8.5 corresponds to 
a frequency width at 150 MHz of about 8 MHz. A bandwidth of 5 MHz would then require an 
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integration time of 500 hours. Such a survey operating over several years would enable a substantial 
range in redshift to be probed. 

5. Conclusions 

The history of the universe between the epoch of recombination and the birth of bright sources 
of radiation by z ^ 5 is currently shrowded in mystery: the absence of radiation sources leaves this 
period a "dark age." Yet spectra of QSOs show not only metal absorption features in their spectra, 
the signature of ongoing and past star formation, but evidence for non-linear density fluctuations 
on large scales as well, a "cosmic web" in the form of the Lya forest. 

We have shown that it may be possible to detect both the development of the cosmic web 
and the epoch of the first radiation sources, whether stars or QSOs, using 21-cm measurements 
from intergalactic gas. Prior to reionization, the radio signal will display emission and/or absorp- 
tion structures both in angle across the sky and in redshift. Fluctuations in 21-cm emission will 
result from non-uniformities in the gas density and from non-uniformities in the distribution of 
Lya radiation sources required to decouple the spin temperature from the CMB temperature. With 
an instrument like the Square Kilometre Array, such fluctuations would be detectable on arminute 
scales, and so reveal whether the reheating is due to galaxies or QSOs. 

Absorption measurements may permit the epoch of the first stellar sources to be determined. 
The continuum radiation from an early generation of stars, redshifted to Lya , will couple the spin 
temperature to the low kinetic temperature of the IGM for a transitory interval before heating the 
gas above the CMB temperature. The result is a detectable absorption signature against the CMB 
that would flag the first epoch of massive star formation in the universe. If reheating were due to 
ultraluminous but sparsely distributed sources (e.g., the first bright QSOs), the resultant patches 
would have a larger scale than any gravitationally-induced inhomogeneities. Because of the very 
large scales involved, there is no difficulty in principle in achieving a sufficiently narrow bandwidth 
and an adequate angular resolution with the SKA or GMRT; the main limitations are the low flux 
levels and the radio background subtraction. Alternatively, a special purpose telescope dedicated 
to scanning the sky in frequency along cold lines of sight through the Galaxy may be used to detect 
the epochs of reheating. Subsequent high resolution measurements with the SKA would provide 
a means of identifying the nature of the first radiation sources and so shedding new light on the 
thermal and reionization history of the universe. 

Support for this work was provided by NASA through ATP grant NAG5-4236 (PM, AM, and 
PT), and by the Royal Society (MJR). We thank the anonymous referee for clarifying the technical 
requirements needed for this project. The results in this paper made use of the Hydra iV-body 
code (Couchman, Thomas & Pearce 1995). 
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A. The Wouthuysen Field mechanism 

Electric dipole radiation (Lya photons) can induce a spin transition through spin-orbit cou- 
pling. The relevant quantum number is the total angular momentum F = I + J, where / is the 
proton spin and J is the total angular momentum of the electron, J = S + L. There are six hyper- 
fine states in total that contribute to the Lya transition. Only four of these, the n = 1 singlet o>Si/2 
and triplet 1S1/2 states (the notation is fLj), and the two triplet n = 2 states 1-P1/2 and 1P3/2 con- 
tribute to the excitation of the 21-cm line by Lya scattering. The selection rule AF = 0, 1 permits 
the transitions S l / 2 -> i-Pi/2; 1P3/2 and i p i/2, 1P3/2 i$i/2, and so effectively qS 1/2 1-5*1/2 
occurs via one of the n = 2 states. Denoting the occupation number of o<Si/2 by n o an d that of 
1S1/2 by ni, the rate equation for no is 

^ = A W (l + jr)n 1 - 3^io^n + A>i - P >o, (Al) 

where the radiation intensity at the 21-cm frequency has been expressed in terms of the antenna 
temperature Tr, I±q = 2A;bTr/Ai . The ratio Tr/T* is the number of 21-cm photons per mode. 
Here, Pq[ and P-fg are the effective excitation and de-excitation rates due to Lya scattering. These 
may be related to the total scattering rate of Lya photons, P a , as shown below. 

It is convenient to label the n = 2 levels as states 2 — 5, from lowest energy to highest. Then 
the external radiation intensity at the frequency Vij corresponding to the i — > j transition may be 
expressed as the antenna temperature , as above. The temperature corresponding to the energy 
difference itself, hvij, is expressed as T-ij = hfij/k-Q. If Aji denotes the spontaneous decay rate for 
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the transition j — > i, then Pq[ and P|g are given by 

p q _ Q T 03 ^30^31 Q To4 ^40^41 , 

^03 ^30 + ^31 J 04 A A0 + All 

and 

pa = ^13 ^30^31 ?14 ^40^41 

10 T 13 A 30 + A31 Ti4 A 40 + An ' 1 ' 

The ratios Aji/A a , where A a is the total Lya spontaneous decay rate, may be solved for using 
a sum rule for the transitions. This states that the sum of all transitions of given nFJ to all the 
n'J' levels (summed over F' and mp>) for a given n'J' is proportional to 2F + 1 (e.g. Bethe & 
Salpeter 1957). There are four sets of downward transitions to n' = 1, J' = 1/2, corresponding to 
the decay intensities I 5 i, ho(= 0), hi, ho, hi, ho, hi, and ho(= 0). These give 

hi _ 5 ho + hi _ ^ ho + ^31 _ g 
ho + hi 3 /30 + hi hi 

Similarly, there are four sets of upward transitions to n' = 2, J' = 1/2 or 3/2, giving the intensities 
^30, ho(= 0), I31, hi, ho(= 0), /40, 4i, and hi- These ratios are 

140 1 730 1 (A5) 



hi + hi 3 I21 + /31 3 

Using hj/h = (9k/gtot){Akj/A a ), where I a is the total Lya decay intensity summed over all the 
hyperfine transitions, and gk = 2F+1 is the statistical weight of level h (gtot is the sum of the weights 
of the upper levels), gives A 2 o/A a = A 50 /A a = 0, A21/A2 = A 51 /A a = 1, A30/A* = A 41 /A a = 
1/3, and Ai/A* = Ao/A* = 2/3. We then obtain P a = 3A a (T^i/T a ) (T a = hv a /k-B), and 
P10 = (4/27)P Q . Statistical equilibrium for Tr = requires Pq\/P\o = n%/no = 3exp(— T*/Tk)- 
Solving the rate equation (|Al]) for Tr > gives equation @, where Tr = Tcmb> and where 
Ts S> T*, Tk 3> T*, and Tr 3> T* have been assumed. 
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Fig. 1. — Predicted rms brightness temperature fluctuations relative to the average at 150 MHz 
for the different cosmologies described in the text, as a function of beam size. A neutral medium 
with spin temperature T$ S> 7cmb has been assumed. The bandwidth is 1 MHz (corresponding to 
a comoving length of 6.5 hr 1 Mpc for Qm = !)■ The data points with error bars show the range of 
values found in a set of iV-body simulations for tCDM (squares) and OCDM (triangles). 
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Fig. 2. — Detectability of IGM fluctuations in 21-cm flux density per beam (in /uJy) at 150 MHz 
{z = 8.5) by a projected facility like the Square Kilometer Array as a function of beam width and 
bandwidth. Dashed lines: Size of the rms fluctuations (0.5^Jy, 2fiJy, and 5/uJy), relative to the 
mean expected in a beam defined by the angular diameter A6 and bandwidth Aza Solid lines: 
Threshold below which the fluctuations are undetectable (S/N < 5) with the SKA for integration 
times of 10, 100, and 500 hours (from top to bottom). Note that this calculation assumes that 
the array is optimally resized (becoming more compact at coarser angular resolution) in order to 
match its resolution to the desired angular scale. 
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Fig. 3. — Left: Radio map of redshifted 21-cm emission against the CMB in a tCDM cosmology 
at z = 8.5. The linear size of the box is 20/i -1 (comoving) Mpc. The point spread function of the 
synthesized beam is assumed to be a spherical top-hat with a width of 2 arcmin. The frequency 
window is 1 MHz around a central frequency of 150 MHz. The color intensity goes from 1 to 6 //Jy 
per beam. For clarity, the contour levels outline regions with signal greater than 4 /iJy per beam. 
Right: Same for OCDM. 
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Fig. 4. — Top: Mean AT = Tj, — Tcmb versus frequency for an angular resolution of 1 arcmin and 
frequency resolution of 1 MHz (solid line), together with the rms fluctuation value (dashed line). 
The IGM is reheated at z t ^ = 9. The strong absorption feature in the mean antenna temperature 
is associated with the fast rise of a Lya continuum background on a time scale « 10 Myr (see text), 
coupling Ts with Tk- Bottom: Corresponding evolution of the kinetic (dashed line), spin (solid 
line), and CMB temperatures (dotted line). 
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Fig. 5. — 21-cm emission and absorption against the CMB from the region surrounding a quasar 
source (lower right corner, in the center of an H n zone), revealed once the IGM is heated above 
the CMB by soft X-rays from the quasar. The linear size of the box is of 20h~ 1 comoving Mpc 
(tCDM), the angular resolution is 2 arcmin, and the bandwidth is 1 MHz. The color levels range 
from —3 /iJy to 3 //Jy per beam. The dark contour demarcates the absorption ring. The quasar 
turns on at z = 8.5 with an ionizing photon luminosity of 10 57 photons s _1 , and is observed after 
7 (left) and 10 (right) Myr. The temperature of the IGM beyond the light radius is assumed to be 
T K = 26mK(l + z) 2 w 2.3 K < Tomb- 



